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Fermi surface models applied to the underdoped cuprates predict the small pocket area to be 
strongly dependent on doping whereas quantum oscillations in YBa2Cu306+ x find precisely the 
opposite to be true — seemingly at odds with the Luttinger volume. We show that such behavior 
can be explained by an incommensurate antinodal Fermi surface nesting-type instability — fur- 
ther explaining the doping-dependent superstructures seen in cuprates using scanning tunneling 
microscopy. We develop a Fermi surface reconstruction scheme involving orthogonal density waves 
in two-dimensions and show that their incommensurate behavior requires momentum-dependent 
coupling. A co-operative modulation of the charge and bond-strength is therefore suggested. 



Identification of the forms of order competing with 
superconductivity and antiferromagnetism in the high- 
T c cuprates remains a considerable experimental chal- 
lenge dJU]. Among possibilities, charge ordering is re- 
ported in several experiments within the underdoped 
regime — namely x-ray diffraction [3], neutron scatter- 
ing [4], scanning tunneling microscopy (STM) [5\ and nu- 
clear quadrupole resonance (NQR) (see Fig.[l}i). Yet 
its extent and relevance are far from understood. It is yet 
to be established whether such ordering participates in 
forming the pseudogap 5 , whether it is inherently uni- 
directional as opposed to bidirectional in nature [7], or 
whether it is caused by a Fermi surface instability [5] as 
opposed to being a biproduct of spin order [3] . 

In the light of recent quantum oscillation [SJ [TU] , elec- 
trical transport [TT] and angle-resolved photoemission 
spectroscopy (ARPES) [12 studies, several Fermi surface 
reconstruction models have been invoked in the under- 
doped cuprates postulating charge (and/or other forms 
of) ordering [TU1 [T3TH3] . A serious problem with all pro- 
posed models, however, is that they predict the pocket 
size to be strongly dependent on the hole doping (e.g. 
dotted and dashed lines in Fig. [l]o) , whereas experiments 
on underdoped YP^CuaOe+z fTt)HH)] find the pocket 
area to change remarkably little over a range of hole dop- 
ings spanning « 3 % [5D] (circles in Fig. [lb). 

In this paper, we show that the near doping- 
independence of the orbit area in underdoped 
YBa 2 Cu30 6 +a; [HI El [20] an d the increasing charge 
modulation period seen with hole doping in STM exper- 
iments on Bi2_j / Pbj,Sr2_ 2 La z Cu06+x [E] can both be 
consistently explained by Fermi surface reconstruction 
resulting from an antinodal Fermi surface nesting-type 
instability (i.e. at [±f,0] and [0,±f] in Fig. [5ji). We 
present a density wave model in which we mimic incom- 
mensurate behavior by considering modulation periods 
A corresponding to different rational multiples of the in 
plane lattice vectors (e.g. A = 7 /2, 4, 13 /3, 9 /2, 5, 6 and 
7). On treating scenarios in which the coupling between 
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FIG. 1: a Charge modulation periods seen using x-rays 
and NQR in YBa2Cu306+^ [6] (large triangles), neu- 
trons in Lao.4Ndo.4Sro.i2Cu04 [3] (pentagon) and STM 
in Bi2-;,Pbj,Sr2- z La 2 Cu06+a; (diamonds), Bi2Sr2CaCu20s+5 
(squares) and Ca2-a ! Na :E Cu02Cl2 (small triangle) taken from 
Ref. [8]. In comparing different materials, we neglect possi- 
ble differences in e(k) [22]. The line and circles show the p 
for each A extracted from the model density-of-states min- 
imum (e.g. Fig. [3J3). b Measured leading YBa2Cu306+a; 
quantum oscillation frequency [16H2U] (circles) compared to 
its strong p-dependence expected in the 4 hole pocket [I] 
(dotted line), Millis and Norman stripe [13] (dot-dash line) 
and fixed A = 4 bidirectional charge |15| (dashed line) mod- 
els, where F = (h/2ire)A c . The present model (solid line) 
uniquely yields a weakly p-dependent F [20] . 



translated bands is uniform (as in a charge-density 
wave [13[ 115] ) or acquires a momentum-dependence (as 
occurs on incorporating a bond-density wave compo- 
nent [2~T]), we find that only the latter leads to a single 
well-defined gap in the electronic density-of-states at 
weak couplings V x ,y <C ^io (where iio is the nearest 
neighbor hopping [22] )■ We show the latter also to be a 
necessary prerequisite for incommensurate behavior, in 
which the electronic structure evolves continuously as a 
function of A. 

We model Fermi surface reconstruction caused by mod- 
ulations of general period A = n /m (in which n and m are 
integers) along the a and/or b lattice directions by diag- 
onalizing a Hamiltonian consisting of nested matrices 
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Here, I„ is an identity matrix of rank n, n! = n for bidirectional order (or n? = 1 for unidirectional order), 
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and £_7'Q x +iQ represents the electronic dispersion 
e(k) [22] subject to translation by multiples of Q x = 
[^,0]andQ„ = [0,^]. 

In the case of a conventional density wave, the normal 
assumption is for the potentials to uniformly couple all 
band crossings subject to a relative translation by 
or Q y such that V x = V x ,o and V y = V Vt o are constants 
in Equation In the case of incommensurate order- 
ing in a two-dimensional lattice, however, the coupling V 
has been found to vary depending on the band crossing 
in question [231 12S] ■ Such behavior is most apparent in 
i?Te3 [25j (owing to its exceptionally large gap), where 
ARPES finds a momentum-dependent V(k) that selec- 
tively couples portions of the Fermi surface subject to 
nesting. 

While the real-space implications of a momentum- 
dependent coupling in the chalcogenides has yet to be 
investigated, in the cuprates it is connected with the pos- 
sibility of bond-strength or bond-current density wave 
ordering [2T] . We find a simple form of the coupling [2~7] , 



nlb r - 



V x (k) = V xfil ^(l-rcosbk y ) 
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in which r adds a bond-strength modulation to an other- 
wise conventional charge-density wave, to prove particu- 
larly effective at reducing the electronic density-of states 
(and consequent free energy) when r sa 1 [2~8]. It does so 
by suppressing V(k) in the regions of the Brillouin zone 
where unnested bands cross [28j . which we demonstrate 
in Fig. [2] by considering the simple case of a unidirec- 
tional modulation Q x — [f^jO] [in which A = 4, n = 4 
and n' = 1 in Equation Q]. 

From Figs. [2)3 and c it is evident that while both uni- 
form (r = 0) and strongly momentum-dependent (r = 1) 
forms of V x open a gap at \k y \ > Sr, where the Fermi sur- 
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FIG. 2: a The unreconstructed Fermi surface at p = 8.5 % 
(black) 22\ together with itself translated by multiples of Q x 
(cyan) and multiples of and Q H (grey) for A = 4. b Recon- 
structed Fermi surface (black) resulting from a unidirectional 
charge modulation (i.e. n' — 1) in which V Xt o ~ 0.3t and r = 
in Equation Q , shown for a quadrant of the extended Bril- 
louin zone, c Same as (b) but with a momentum-dependent 
V x in which we choose r = 1 [28]. d The calculated electronic 
density-of-states (DOS) for the unreconstructed band [22j ex- 
hibiting a van-Hove singularity near e w — 2t. e The calcu- 
lated DOS (black line) for r = in (b). f The calculated 
DOS (black line) for r — 1 in (c). Red lines in (e) and (f) 
are the corresponding DOS calculated for concurrent charge 
modulations along a and 6 (i.e. such that n' — n) in which 
we assume V x ,o = Vy,o (by no means a required constraint). 



faces are nested by Q^, the latter does so without split- 
ting the open Fermi surfaces at k y w ±|r. The splitting 
in Fig. [2)3 occurs concomitantly with an additional gap 
in the electronic density-of-states at e « -1.8 t in Fig. [2^ 
and a slightly weaker ordering gap at the Fermi energy 
(e F sa — t) than in Fig. [2]f. A large V x at \k y \ sa |j is 
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therefore energetically unfavorable . The momentum- 
dependent V x (i.e. r m 1) avoids unfavorable splittings 
and gaps, leaving the remaining open Fermi surfaces at 
k y « ±^ amenable to a secondary Fermi surface insta- 
bility of wavevector Q y = [0, which can further lower 
the density-of-states (and consequently the electronic en- 
ergy) by forming a concurrent modulation along b (red 
line in Fig. [2)f) [where n' = n = 4 in Equation (JlJ in the 
case of bidirectional ordering]. By contrast, the splittings 
caused by a uniform potential (i.e. r = 0) mutually dis- 
rupt nesting for both Q x and Q y in the case of bidirec- 
tional ordering, leading to an energetically unfavorable 
higher density-of-states consisting of multiple peaks and 
valleys in the vicinity of the Fermi energy (red line in 
Fig. §3). 

On extending the bidirectional ordering density-of- 
states calculation to different periods in Fig. [3j we con- 
tinue to find a well defined single gap with a broad deep 
minimum only in the case of momentum-dependent cou- 
pling (see Fig. J3Jd) , pointing to its continuous evolution 
with A. In the case of a uniform coupling (see Fig. [3£i), 
by contrast, the multiple peaks and valleys vary discon- 
tinuously with A. 





FIG. 3: a Electronic density-of-states (DOS) calculated for 
bidirectional charge modulations with periodicities corre- 
sponding to different multiples A of the a and b lattice vectors 
as indicated, assuming uniform couplings V x ,o = V Vl o — 0.3t 
in which r = 0. b Same as (b) but assuming momentum- 
dependent couplings in which r = 1 in Equation (pi. Curves 
have been offset for clarity. The dashed line in (b) indicates 
the minimum in the DOS near to which the Fermi energy is 
likely to be located. 



Thus, by generating a deep wide gap in the density-of- 
states whose form and location in energy shifts continu- 
ously with A, momentum-dependent coupling provides an 
incentive for incommensurate behavior in which A adjusts 
itself in a continuous fashion so as to lower the electronic 
energy. Because the electronic energy in an itinerant pic- 
ture is minimized by having the Fermi energy situated 
within a broad deep gap in the density-of-states, the A- 
dependent gap provides an explanation for the evolution 
of the periodic structures seen in STM experiments as 



a function of doping [8] . The location of the minimum 
(identified by the dot-dashed line in Fig. [3)3) enables us to 
estimate the hole doping p at which each period is most 
likely to be stable (plotted in Fig. []Ji). Using these dop- 
ings and assuming Luttinger's theorem we calculate 
the corresponding Fermi surfaces in Fig. [4j whose forms 
consist of a single electron orbit (located close to the 
nodes) consistent with experimental observations [151130]. 
Momentum-dependent coupling enables such a pocket to 
exist for weaker couplings than in ref. [15] and to per- 
sist essentially unchanged as a function of doping. Most 
importantly, the near p-independent area (solid line in 
Fig. [TJd) reproduces experimental observations (circles) . 
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FIG. 4: a, b, c and d Reconstructed Fermi surface for se- 
lected A's in Fig. (3ja when the Fermi energy is situated at 
the minimum in the density-of-states, with the corresponding 
hole doping given. Solid red lines indicate the fc-space area 
of the Aa x Ab superstructure, while dashed red lines indicate 
the n 2 -fold reduced Brillouin zone [which coincides with the 
superstructure in (a) and (d)]. In (b) and (c), a magenta line 
is used to trace the path of the electron-like orbit that occurs 
in strong magnetic fields, ho ~ 100 meV |22| produces an 
effective mass consistent with experiments. 



The sub-gaps occurring at the intersections of the 
electron orbits in Figs. [4}d and c are small enough 
[A^ uh /BF(he/m*) 2 < 1 provided V XtV < ho] to be com- 
pletely broken down [31] in magnetic fields of the strength 
required to see magnetic quantum oscillations [5RTT1 UM 
[T71 [3"0] — giving rise to a single orbit (thick magenta 
line) in strong magnetic fields. The sub-gaps neverthe- 
less imply the absence of a simple (A-independent) rela- 
tionship between the quantum oscillation frequency F c = 
(H/2ne)A e and the frequency F L = |Fbz correspond- 
ing to the Luttinger hole doping (where Fez = h/eab 
is the unreconstructed Brillouin zone frequency). Only 
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when the density-wave is 'accidentally' commensurate 
such that sub-gaps do not occur (e.g. A = 4 or 5) can 
adherence to Luttinger's theorem [29] be easily verified 
in quantum oscillation experiments. In Fig. [4^i, for ex- 
ample, Fl = F\ — F c (where F\ = Fbz/A 2 is the Xa x Xb 
superstructure frequencsy), while in Fig. [4jd it is given 
by F L = lF x - F c . 

Finally, we turn to aspects of momentum-selective den- 
sity waves that may potentially be reconciled with the 
unidirectional behavior of charge ordering noted in the 
cuprates [7j. While closed Fermi surface pockets require 
charge modulations to occur concurrently along a and b 
lattice directions (in the absence of other orders [E]), the 
superposition of their ordering gaps in Fig. |2jF (red line) 
implies the absence of a significant energy penalty (or in- 
teraction) associated with their coexistence — in contrast 
to the uniform case in Fig. [2^ where such a superposi- 
tion does not occur. Given the implied independence of 
the modulations along a and b, underlying anisotropics in 
the electronic structure (such as that caused by the pres- 
ence of chains in YE^CuaOe+z [23 ) will likely produce 
anisotropies in V XtV fi, A, r and the onset temperature. In 
the present simulations, we find a Fermi surface consist- 
ing solely of an electron pocket to remain robust against 
an anisotropy V^o/K/.o as large as 4. 

In conclusion, we present a model that explains the 
lack of a detectable doping-dependence of the quantum 
oscillation frequency in underdoped YE^CusOe+x (i.e. 
Fig. [TJd [20]). By considering rational values of A, we 
develop what is in essence an incommensurate model for 
co-operative charge- and bond-density wave ordering in 
the cuprates [U] [31] — here driven by a Fermi surface 
instability at the antinodes. By incorporating a (possi- 
bly dominant 28 ) bond-density wave component [2"T] . 
the size of the periodic potential required to produce 
a single pocket with a small residual density-of-states 
is greatly reduced (i.e. V x , y ,o ^ 0.05iio) [IE] relative 
to other models [T3TU5| . A key strength of the present 
model is its ability to reconcile doping-dependent quan- 
tum oscillations [T6H20"] [30] with the doping-dependent A 
seen m STM and other experiments [3 [6] [8] (i.e. Fig.[l]b), 
the negative Hall and Seebeck coefficients over a broad 
range of dopings seen in transport [11] and particle- 
hole symmetry breaking reported at the antinodes in 
ARPES [2] [12] — all while maintaining compliance with 
Luttinger's theorem [29] . 
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